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ABSTRACT 

Chandra X-ray observations of the nearby brightest cluster galaxy M87 resolve the hot gas 
structure across the Bondi accretion radius of the central supermassive black hole, a mea¬ 
surement possible in only a handful of systems but complicated by the bright nucleus and 
jet emission. By stacking only short frame-time observations to limit pileup, and after sub¬ 
tracting the nuclear PSF, we analysed the X-ray gas properties within the Bondi radius at 
0.12 — 0.22 kpc (1.5 — 2.8 arcsec), depending on the black hole mass. Within 2 kpc radius, 
we detect two significant temperature components, which are consistent with constant val¬ 
ues of 2 keV and 0.9 keV down to 0.15 kpc radius. No evidence was found for the expected 
temperature increase within ^ 0.25 kpc due to the influence of the SMBH. Within the Bondi 
radius, the density profile is consistent with pocr^^ The lack of a temperature increase inside 
the Bondi radius suggests that the hot gas structure is not dictated by the SMBH’s potential 
and, together with the shallow density profile, shows that the classical Bondi rate may not 
reflect the accretion rate onto the SMBH. If this density profile extends in towards the SMBH, 
the mass accretion rate onto the SMBH could be at least two orders of magnitude less than 
the Bondi rate, which agrees with Faraday rotation measurements for M87. We discuss the 
evidence for outflow from the hot gas and the cold gas disk and for cold feedback, where gas 
cooling rapidly from the hot atmosphere could feed the cirumnuclear disk and fuel the SMBH. 
At 0.2 kpc radius, the cooler X-ray temperature component represents ~ 20% of the total X- 
ray gas mass and, by losing angular momentum to the hot gas component, could provide a 
fuel source of cold clouds within the Bondi radius. 
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1 INTRODUCTION 


Accretion onto the supermassive black hole (SMBH) at the centre 
of the giant elliptical galaxy M87 is powering relativistic jet activity 
that hea ts the surrounding hot clus ter atmosphere and stifles galaxy 
growth ^Forman et alj|200^l2007l) . Deep Chandra X-ray observa¬ 
tions of the intracluster medium have found large cavities filled 
with radio emission where the X-ray gas ha s been displaced by the 
expansion of the je t l Bohringer^DfLl 1.9951: Bicknell & BegelmanI 
Il996l : lYoung et alj|2002l: Iporman et alJbOOSL bOolT This energy 
input is likely suppressing the rapid cooling of the dense, low tem¬ 
perature X-ray gas at the centre of the cluster, which would oth¬ 
erwise flood the brightest cluster galax y (BCG) with cold molec - 
ular gas and ensuing star formation (eg. IPeterson & Fabiaj|200d) . 
In general, this feedback from the central AGN is thought to be 
the essential mechanism in models of galaxy formation that trun- 
cates galaxy growth by restricting gas cooling and star fo rmation 
dCroton et al.ll2Q0^ : iBower et alj2006l : iHopkins et alj|2003) . 


* E-mail: hrr27@ast.cam.ac.uk 


The central AGN is located in the densest region of the clus¬ 
ter’s hot atmosphere and therefore classical Bondi accretion of the 
hot gas is of t en considered as a fuelling mechanism for feedback 
(lBondilll95l : lBaganoff et alj2003l : lDi Matteo et alj2003h . The gas 
will be accreted by the SMBH if it lies within the Bondi radius, 
where the gravitational potential of the SMBH dominates the ther¬ 
mal energy of the hot atmosphere. The Bondi radius is given by 
r'B = 2 GMbh/cj, where Mbh is the mass of the black hole and Cj 
is the sound speed of t he hot gas. Using Chandra observations of 
M87, iDi Matteo et alj j2003h found that the Bondi accretion rate 
is able to provide a sufficient supply of fuel to power the AGN. 
However, Bondi accretion of the hot gas m ay be insufficient to 
fuel the most powerful radio jet outbursts (eg. lRaffertv et al.ll20Q6l : 
iMcNamara et alJ I2OI ih . Cold gas accretion from circumnuclear 
disks of atomic an d molecular gas may be r equired to supplement 
the fuel supply feg. lPizzolato & Sokeilf200^ . 

For accretion at the Bondi rate, the radiative output from 
the nucleus in M87 is orders of mag nitude below the expecta¬ 
tion from standard th i n disk accretion ( [Fabim & CanizarM i 19881 : 
iFabian & Reed 1 19951 : IdI Matteo et alj 120001) . The black hole at 
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Table 1. Details of the Chandra observations used for this analysis and best-fit nuclear spectral model parameters. 


Obs. ID 

Date 

Exposure 

Frame time 

«H 

r 

Flux (2-lOkeV) 



(ks) 

(s) 

( 10^2 cm- 2 ) 


(10^*^ erg ) 

352 

29/07/2000 

29.4 

3.2 

- 

- 

- 

1808 

30/07/2000 

12.8 

0.4 

0.08±0.01 

2.37 ±0.07 

0.70 ±0.04 

11513 

13/04/2010 

4.7 

0.4 

0.06 ±0.02 

2.35 ±0.08 

^•^^-0.09 

11514 

15/04/2010 

4.5 

0.4 

0.05 ±0.02 

2.15 ±0.09 

1.4±0.1 

11515 

17/04/2010 

4.7 

0.4 

0.04 ±0.02 

2.21 ±0.08 

1.6±0.1 

11516 

20/04/2010 

4.7 

0.4 

0.02 ±0.02 

2.02 ±0.08 

1.7±0.1 

11517 

05/05/2010 

4.7 

0.4 

0.08 ±0.02 

2.34 ±0.08 

1 54+010 
^•^^-0.09 

11518 

09/05/2010 

4.0 

0.4 

0.06 ±0.02 

2.3±0.1 

1 17 + 0.10 
^•^'-0.09 

11519 

11/05/2010 

4.7 

0.4 

0.08 ±0.02 

2.4±0.1 

1.02 ± 0.08 

11520 

14/05/2010 

4.6 

0.4 

0.11 ±0.02 

2.5±0.1 

^•^^-0.07 

13964 

05/12/2011 

4.5 

0.4 

0.07 ±0.02 

9 T 1 + 0.10 

^•■^^-0.09 

1 17+0.09 
^•^'-0.08 

13965 

25/02/2012 

4.6 

0.4 

0.04 ±0.02 

2.11 ±0.09 

1.5±0.1 

14973 

12/03/2013 

4.4 

0.4 

0.05 ±0.02 

2.3±0.1 

1.18±0.09 

14974 

12/12/2012 

4.6 

0.4 

0.05 ±0.03 

2.2 ±0.1 

1.17±0.09 

Stack 


54.7 

0.4 

0.057 ±0.006 

2.25 ±0.03 

1.30 ±0.03 


the centre of the Milky Way, Sgr A*, is a good example of a 
low luminosity AGN with a radiative output several orders of 
magnitude below that implied by the Bon di accretion rate and 
the standard 10% radiative efficiency (eg. iBasanoff et al.l 1200 3l : 
lOenzel et alj l2010h . This observed radiative inefficiency of the 
accretion disk is consistent with an advection-dominated accre¬ 
tion flow where the gas is unable to radiate the gravitational po¬ 
tential e nergy locally an d the trapped ener g y is instead advected 
inward jlchimarul Il977l: iRees et al ] 1 19821: iNarayan & yH 1 1994 
iNarayan & McClintockl2008l) . Alternatiyely. the mass inflow at the 
Bondi radius may not reach the event horizon of the SMBH but is 
instead ejected from the black hole ’s sphere of influence (ADIOS; 
[Blandford & Bege lmanlll999| . 2004 or circulates in convectiv e ed- 


dies (CDAF: lNaravan et al.ll2000l : Ouataert & Gruzin^l2000h . 

It has proven difficult to distinguish between these possibili¬ 
ties, or indeed to determine the fuelling mechanism of the AGN, 
because it is only possible to resolve the hot gas inside the Bondi 
radius in a handful of systems. Deep Chandra observations of 
two of the best targets , Sgr A* and NGC3115 jBaganoff et all 
l2003l:IWong etalj|201lh . are consistent with flat den sity profiles 
with inside the Bondi radius dWang et al.l2013l : IWong et'd] 

I20l4 . This is consistent with recent numerical simulations which 
show outflows in the form of winds off the accretion flow preyent- 
ing all but a fraction of the Bondi r ate from reaching the SMBH 
(eg. lYuan et aDl2012l : iLi et al.ll2013h . However, whilst such a low 
mass accretion rate is adequate for these quiescent systems, it is 
insufficient to fuel powerful radio-jet outbursts without requiring 
effici encies much greater than 100% dNemmen & TchekhovskovI 
|2014 . Such high efficiencies are theo retically achievable if the 
black hole spin energy can be extracted l Tchekhovskoy_eL^ .1201 ll: 
Tchekhoyskoy & McKinney! l2012l : iMcKinney et al.ll2012h but 
tematiyely, the inflow rate may be significantly increased by 
cold gas accretion, which compensates for the outflowing mass 
dPizzolato & Sokeill200.5l : iGasnari et alJl2013ll . 

At a distance of only ~ 16 Mpc, M87 at the centre of the 
Virgo cluster is one of few systems for which Chandra can re- 
solye the Bondi radius at ^ 2 arcsec (~ 0.15 knc: IPi Matteo et ^ 
I2OO3I : iGarcia e^ ^l2010h . H osting a powerfu l radio jet outburst 
and a circumnuclear g as disk dpord et alJll994 iHarms et alJll994 
iMacchetto et alJl997tl . M87 provides an excellent comparison sys¬ 
tem to NGC3115 and Sgr A* to determine the properties of the 


putative accretion flow. However, the bright nuclear point source 
and jet knots make it challenging to extract the properties of the 
underlying cluster gas across the Bondi radius. Early Chandra ob¬ 
servations of M87 in 2000 showed that the bright nucleus and jet 
emission would not be piled up in a 0.4 s frame t ime observation 
I Wilson & Yangll2002l : IPerlman & Wilsor]l2005h . IPi Matteo et al.l 
1 2 OO 3 I) used this short 12.8 ks observation to determine the gas 
properties using single temperature spectral models in to 2 arcsec 
radius. In 2005, major flaring of the jet knot HST-1, located only 
0.85 arcsec from the nucleus, th en resulted in stro ng pileup even 
in short frame-time observations dHarris et alj200d) and prevented 
follow up studies. However, HS T-1 has slow l y decl ined in bright¬ 
ness since the peak in 2005 and [Harris et al.l 1201 ih found that by 
2010 the intensity had dropped to the level observed in 2000. 

Here we stack twelve 5 ks monitoring observations of M87 
taken since April 2010 (PI Harris) and subtract the bright nuclear 
emission to determine the density profile inside the Bondi radius 
and study the temperature structure on these scales. Section de¬ 
scribes the observations, data reduction and subtraction of the nu¬ 
clear emission. In section[3] we determine the temperature and den¬ 
sity structure within the Bondi radius and calculate the cooling and 
dynamical time profiles. We discuss accretion of the hot gas atmo¬ 
sphere by the SMBH and the evidence for hot and cold gas out¬ 
flows in section |4] before presenting our c onclusions in sect ion [5] 
We assume a distance to M87 of 16.1 Mpc dTonrv et alj200lh to be 
consistent with earlier analyses of the Chandra datasets. This gives 
a linear scale of 0.078 kpc arcsec^ * and ~ 2 — 5 x 10^ gravitational 
radii per arcsec, depending on the black hole mass (section[T4j. All 
errors are let unless otherwise noted. 


2 CHANDRA DATA ANALYSIS 

The centre of M87 has been observed regularly with Chandra since 
2002 to monitor jet activity, understand the emission mechanism 
and determi ne the site of the TeV flaring dHarris et alJl2003ll2006L 
I 2 OO 9 I I 2 OI ih . For this study of the central cluster gas properties, 
we have selected all archival Chandra observations of M87 taken 
since the intensity of HST-1 dropped back to its 2000 level and 
after the TeV flaring ev ent in early April 2010 dOng & Mariottil 
l201(t lllarris et al]|201 ih . These ACIS-S datasets were also com- 
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Figure 1. Left: Background-subtracted surface brightness profiles in the energy range 0.5 — 7.0 keV for obs. ID 1808 and the summed obs. IDs 11513 to 
14974. The ChaRT simulations of the nuclear flux for each observation are shown by the open symbols. Right: Same as left but with the nuclear flux subtracted 
from each profile. The radial range for the Bondi radius is shown by the vertical dashed lines. 


pared with two of the ACIS-S observations taken in 2000 to test the 
subtraction of the nuclear emission (obs. ID 1808) and determine 
the gas properties to large radii for the deprojection analysis (obs. 
ID 352). We have not included the short, < 1 ks, exploratory ob¬ 
servation, which was taken to determine the nuclear count rate and 
therefore the frame time suitable for the later observations (obs. ID 
351). 

All datasets were analysed with CIAO version 4.6 and 
CALDB version 4.6.2 supplied by the Chandra X-ray Center 
jFruscione et alj|2006ll . Level 1 event files were reprocessed with 
the latest gain and charge transfer inefficiency correction applied 
and then filtered to remove cosmic ray events. Background light 
curves were extracted from regions free of point sources and used 
to remove periods affected by flares from the event files. The com¬ 
bined dataset was also examined for any periods with high back¬ 
ground count rates. No flares were found in the short frame-time 
observations and only two short flares occurred during obs. ID 352. 
The final cleaned exposure times are given in Table[T] 

The cluster emission from M87 extends across the entire 
Chandra field of view, therefore blank-sky background observa¬ 
tions were used to subtract the background from images and spec¬ 
tra. The appropriate blank-sky background dataset was processed 
identically to the event file, reprojected to the same sky position and 
normalized so that the count rate matched that of the event file for 
the 9.5 — 12 keV energy band. Point sources were identified using 
a hard band (3 — 5 keV) image for each observation and excluded 
from the analysis. The shallow datasets from 2010 to 2013 were 
also reprojected to a common position and combined to ensure that 
no faint point sources were misse d. This method iden tified and ex¬ 
cluded all point sources found bv ijordan et alj (l2004ll . The remain¬ 
ing contaminating flux from unresolved low mass X-ray binaries 
and the stellar population (cataclysmic variables and coronally ac¬ 
tive binaries) is insignificant at only ~ 1% of the total emiss ion in 
M87 in the 0.5 — 2 keV energy band teevnivtsev et al.ll20Q8h . The 


cluster emission dominates at energies below 10 keV and the point 
source emission only becomes significant at higher energies. The 
emission from the jet was also carefully excluded. 

2.1 Subtracting the nuclear PSF 

Within 1 arcsec, the nuclear point source is more than an order of 
magnitude brighter than the cluster emission. Although the clus¬ 
ter emission dominates beyond ~ 3 arcsec, the Bondi sphere is lo¬ 
cated within this region. It was therefore necessary to accurately 
model and subtract the nuclear point spread function (PSF) to de¬ 
termine the gas properties across the Bondi radius. This analy¬ 
sis is described in detail in the appendix, including an evaluation 
of the low-level pileup, the nuclear spe ctral model and sim ula¬ 
tions of the PSF. In summa ry, following iRussell et alj fcOlOh and 
ISiemiginowska et"^ i2010t) . we used the best-fit power law model 
for the nuclear spectrum (Table [T]( and the Chandra ray-tracer 
(ChaRT: Icarter et^l2003h to simulate the nuclear PSF (photons 
scattered by the Chandra mirrors). This simulation was then pro¬ 
jected onto the ACIS detector with MARX and used to subtract the 
PSF contamination from the cluster emission at each radius within 
the cluster core. We assume that all emission within 1 arcsec is 
from the AGN. However the gas emission is likely to be 5% of 
the total, depending on the slope of the gas density profile (section 
A2). 

The simulation was tested by comparing the PSF-subtracted 
surface brightness profiles for two separate epochs (obs. ID 1808 
and the stacked 11513 to 14974 observations). Fig. [T] (left) shows 
the total emission profiles and the predicted contribution of the PSF, 
which has varied significantly in flux between the two epochs. Fig. 
[T] (right) shows the subtracted profiles, which incorporate the un¬ 
certainty on the nuclear flux (Table [T]l. The flux of the PSF sim¬ 
ulation was increased to compensate for the reduction due to mild 
pileup within 1 arcsec, which causes an oversubtraction of the AGN 
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emission in this region (section A3). The subtracted profiles agree 
within the errors showing that the majority of the variable nuclear 
emission has been removed and only the cluster emission remains. 
Although the ChaRT/MARX simulation is not expe cted to provide 
an ex act prediction of the PSF to large radii (eg. iGaetz & Jeriu j 
l2005h . the nuclear flux is only significant to ~ 3 arcsec radius and 
the excellent agreement between the subtracted profiles in Fig. [T] 
(right) suggests this method is sufficiently accurate for our analy¬ 
sis. 

The combined effect of the minor inaccuracies in the nu¬ 
clear PSF subtraction due to the cluster contamination, mild pileup, 
nuclear spectral model and effective area reduction (section A3) 
was determined by comparing the subtracted profiles from these 
different epochs. In Fig. [T] the cluster surface brightness in the 
1—2 arcsec radial bin, where the uncertainty due to the PSF sub¬ 
traction is largest, agrees within the errors. This suggests that the 
difference in flux between the two epochs is less than ~ 15%. The 
nuclear PSF flux is significantly lower in obs. ID 1808 therefore it 
is more likely that the PSF contribution has been underestimated in 
the summed observations. To be conservative, we have estimated 
the impact on the measured gas properties if the PSF flux is sys¬ 
tematically higher by 20%. Variation in the photon index by ±0.2 
and absorption by ±0.05 x 10^^ cm^^ produced uncertainties in 
the measured temperature and density values less than the statisti¬ 
cal error and therefore is not significant in comparison to the uncer¬ 
tainty in the nuclear flux. 


2.2 Spectral analysis and deprojection 


Fig. |2] (left) shows the complex cavity structures and bright blobs 
and filaments of cool gas located at the ce ntre of M87 . Thes e 


gas structures have b e en studied in detail b y [)foung_eLaL 
iForman et a D (l2005h. iForman et al (l2007h . iMillion e~ 


Werner et al.l ( iTOlOli . Clear depressions are observed in the X- 


and 

ray emission to the west of the nucleus beyond the jet, immedi¬ 
ately south-east of the nucleus and to the south and east of that 
(~ 20 arcsec from the nucleus). These depressions are cavities 


carved out by the radio 
rounding atmosphere (e 

jet and buoyantly risin 

g through the sur- 

g. Churazov et al.l 200C 

; McNamara et al.l 

iTOOdlFabian et ahlTOOCl: 

McNamara & Nulserj2007h. Fig.l2l(right) 


shows the angular variation in surface brightness around the nu¬ 
cleus. The inner cavity to the south-east can be clearly seen in 
the 140 — 234° sector. This cavity is surrounded by a rim of cool 
gas that is visible as a bright ridge of emission. The X-ray sur¬ 
face brightness to the north and south of the nucleus is clumpy 
and several point sources are visible, which were excluded from 
the profiles. However, there is less angular variation in the surface 
brightness inside ~ 4 arcsec radius. We have therefore proceeded 
with azimuthally-averaged density and temperature profiles. 

Spectra were extracted from circular annuli extending to a 
radius of 14 arcsec centred on the nucleus for obs. IDs 11513 to 
14974. The radial range from these observations was limited by the 
subarray selection and therefore spectra were also extracted from 
obs. ID 352 to extend the radial profiles to ~ 300 arcsec. Although 
there has been a significant decrease in the effective area at low 
ene rgies between these observations due to the contaminant build 
up dMarshall et al.l200^ . the surface brightness profiles are consis¬ 
tent in the region of overlap (section A3). Background spectra were 
generated from the blank sky backgrounds for each source spec¬ 
trum and appropriate responses and ancillary responses were pro¬ 
duced. The spectra were restricted to the energy range 0.5 — 7 keV 
and grouped with a minimum of 25 counts per energy bin. For obs. 


IDs 11513 to 14974, the source was located at a similar position 
on the same chip and therefore the responses were comparable. We 
summed together the spectra and background spectra for each an¬ 
nulus and averaged the response files with a weighting determined 
by the fraction of the total counts in each observation. The summed 
spectra contained a minimum of 3000 counts in each of the nar¬ 
rowest 2 arcsec radial bins and over 40,000 in the broader radial 
bins. 


The summed spectra covering radii out to 14 arcsec and 
the obs. ID 352 spectra from 14 to 300 arcsec radius were de- 
projected with t he model-independent spectral deproiection rou¬ 
tine DSDEPROJ jSanders & Fabianll2007l : iRussell et ^l2008h . As¬ 
suming spherical symmetry, DSDEPROJ takes the background- 
subtracted spectra and performs a geometric deprojection starting 
from the outermost annulus and subtracting the projected emissio n 
off each successive annulus dFabian et alJll981uKriss et alJll983l) . 
The routine was modified to include a correction for the regions 
excluded due to point sources and a correction for the difference in 
effective area between obs. ID 352 and the summed spectra. 

The resulting deprojected spectra were each fitted in 
XSPEC with absorbed single and two-temperature APEC m odels 
jBalucinska-Church & McCammon|[l9^ : ISmith et al.ll200i1) . The 

cluster redshift was fixed to z = 0.0044 and the absorption of 
the c luster gas componen t was fixed to the Galactic column den- 
sitv dKal berla et aL||2005l), consistent w i th previous analyse s (eg. 
IForman et al. 20071 : Million et aD l20ld: l^mer et alJl201^ . The 
temperature, metallicity and normalization parameters were left 
free and the ;{;^-statistic was minimised for the spectral fit. For 
the two-temperature spectral fits, the metallicity was tied between 
the two components. An additional absorbed power-law compo¬ 
nent was included to account for the nuclear PSF contribution 
in the spectral fit to the summed spectra. If the power-law pa¬ 
rameters are left free, the contribution of the nuclear PSF is in¬ 
correctly interpreted in the spectral fitting as hot cluster emis¬ 
sion. The parameters of the power-law component were therefore 
fixed to the best-fit values determined from s pectral fitting to de- 
projected ChaRT/MARX si mulated spectra (eg. lRussell et ahlEzOld : 
ISiemiginowska et'^l2010h . We also generated a deprojected elec¬ 
tron density profile with finer radial binning from the PSF and back¬ 
ground subtracted surface brightness pr ofile incorporating the tem¬ 
perature and metallicity variations (eg. lCavagnolo et ^l2009h . Fi¬ 
nally, the use of the summed spectra was verified by comparing the 
spectral fit to the summed projected spectra with a simultaneous 
fit to all of the individual projected spectra in XSPEC. The best-fit 
parameters for each case were found to agree to within the errors. 


3 RESULTS 

Using a ChaRT/MARX simulation to subtract the nuclear PSF emis¬ 
sion, we are able to determine the gas density down to 1 arcsec 
radius from the nucleus and the temperature down to 2 arcsec ra¬ 
dius. Fig. 13 shows the deprojected single temperature best-fit re¬ 
sults using radial bins of 2 — 4 arcsec width within the central 
14 arcsec (1 kpc) region. The deprojected temperature drops from 
1.88+o;o^ keV at 5 kpc (64 arcsec) to 0.91+o;ii keV at 0.2 kpc 
(2.5 arcsec ). The coolest X-ra y gas is concentrated in the cluster 
centre (eg. lWerner et alfeoiol) producing a particularly steep drop 
in temperature at 1.6 kpc (20 arcsec). The bright, cool rim of the SE 
cavity is located at a radius of 1 kpc and coincides with an increase 
in the gas density by a factor of 2. The density drops at a radius 
of ~ 0.5 kpc due to the SE cavity and this produces larger errors 
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Figure 2. Left: Exposure-corrected image showing the sectors used to compare the angular variation in the cluster surface brightness around the nucleus. The 
excluded jet region is shown by the dashed black box. Right: Background-subtracted surface brightness profiles in the energy band 0.5 — 7.0 keV extracted 
from these sectors. There is a cavity in the sector from 140 — 234“ which complicates the deprojection analysis. The radial range for the Bondi radius is shown 
by the vertical dashed lines. 


on the deprojected profiles due to the reduced number of counts 
in this region. The apparent increase in the deprojected tempera¬ 
ture at 0.3 — 0.8 kpc radius is therefore not significant and the pro¬ 
file is consistent with a modest decline in temperature towards the 
cluster centre. However, it is also clear that a single temperature 
fit is not adequate and multi-temperature gas is present requiring 
a two-temperature spectral model. The metal abundance peaks at 
1.1 ± 0.1 Zq at a radius of 4 kpc and then appears to drop rapidly in 
the cluster centre. This drop can be caused by the Fe bias that occurs 
when a sing le temperatu re model is fitted to a multi-temperature 
medium (eg. [Buotell200(t) . 


3.1 Two temperature models 

Fig. El shows the deprojected temperature and metallicity profiles 
produced by fitting two-temperature models to the spectra extracted 
in broader radial bins. Within 2 kpc, a second temperature com¬ 
ponent provides a better representation of the data with an F-test 
probability < 0.01% that the second component is detected by 
chance. The best-fit temperature values are 1.98 ±0.09 keV and 
1.00 ±0.02 keV. The metallicity now appears to be consistent with 

1 Zq in the cluster centre. By fixing the metallicity to 1 Zq, a two- 
temperature model was also fitted to the spectra extracted from the 
cluster centre, within 2 kpc, in narrower radial bins. Fig. |4] shows 
the resulting deprojected profiles where multi-temperature gas was 
detected significantly within 2 kpc. In addition to the metallicity 
constraint, the temperature parameters were fixed to 0.9 keV and 

2 keV in the 0.3 — 0.45 kpc radial bin where the cavity increased 
the uncertainties. Only the APEC normalization parameters were 
left free for this region. The gas temperatures for the two compo¬ 
nents appear roughly constant within 1 kpc, although there could 
be a modest decrease in the lower temperature component. No evi¬ 
dence was found for a temperature increase at small radii due to the 


gravitational influence of the black hole, although this may occur 
at smaller radii within 0.15 kpc. 

The temperature profile down to a radius of 2 arcsec was not 
significantly affected by the estimated uncertainty in the PSF sub¬ 
traction. In sectionl^ the PSF-subtracted cluster surface brightness 
profiles from obs. ID 1808 and the summed obs. IDs 11513 to 
14974 were compared and found to agree within the errors. How¬ 
ever, the cluster surface brightness in the 1—2 arcsec radial bin is 
also consistent with a 15% difference between the two epochs, 
which may be due to minor inaccuracies in the PSF subtraction 
such as an underestimation of the pileup fraction. The nuclear PSF 
flux is significantly lower in obs. ID 1808 therefore it is more likely 
that the PSF contribution has been underestimated in the summed 
observations. To be conservative, we have therefore estimated the 
impact on the measured gas properties if the PSF flux is higher by 
20%. This will increase the oversubtraction within 1 arcsec where 
the flux has been reduced by pileup but could improve the sub¬ 
traction of the PSF wings. No significant change in the temperature 
profile was found for this increase in the PSF flux. We also note that 
no significant variation in the temperature was found if the PSF flux 
was instead reduced by 20%. Within 2 arcsec, the nuclear PSF dom¬ 
inates over the cluster emission and the uncertainty in the faint hard 
X-ray emission above 3 keV produced a large uncertainty in the gas 
temperature. The gas temperature profile is therefore restricted to 
radii greater than 2 arcsec. 


3.2 Density profile 

The coolest component of the multi-temperature gas in M87 likely 
consists of small blobs and filaments whilst the 2 keV compo¬ 
nent fills the volume, excluding the cavities (I Young et al.ll200j : 
ISnarks et'n]|2004 IWemer et al.]|2010h . By assuming that the dif¬ 
ferent gas phases are in pressure equilibrium, the filling factor of 
the cooler gas component can be calculated and used to deter- 
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Figure 3. Deprojected temperature, metallicity and density profiles com¬ 
paring a single temperature spectral model (solid circles) and a two temper¬ 
ature spectral model (open crosses). The radial range for the Bondi radius 
is shown by the vertical dashed lines. 


mine the gas density. Note that we have ignored magnetic pressure, 
which is likely requ ired to support extended filaments of cold gas 
dFahian et al.ll2008h . The gas density of both temperature compo¬ 
nents continues to rise towards the cluster centre and the profile 
appears to steepen inside 0.5 kpc radius. A similar increase is seen 
in Fig.[^for the single temperature density profile. To determine if 
this steepening is significant, we have also produced a deprojected 
density profile with narrower 1 arcsec radial bins from the PSF and 
background-subtracted surface brightness profile (see section l2^ . 

Fig-[5]shows the higher spatial resolution density profile which 
extends in to 1 arcsec. The density profile steepens within 0.3 kpc 
and peaks at a density of 0.62 ± 0.05 cm^^ . Measurement of the 
gradient is however complicated by the cavity at ~ 0.5 kpc radius, 
where the emitting volume is overestimated and therefore the den¬ 
sity is underestimated. Due to the cavity, the majority of the emis¬ 
sion at this radius is projected. Therefore the density value is also 
sensitive to the PSF subtraction despite this accounting for only 
~ 2% of the total emission. The nuclear PSF accounts for roughly 
half of the total emission in the 1 — 2 arcsec radial bin and therefore 
any inaccuracy in the subtraction could have a significant effect 



10 100 
Radius (arcsec) 


Figure 4. Deprojected temperature, metallicity and density profiles for a 
two temperature spectral model. The finer radial binning was achieved by 
fixing the metallicity inside a radius of 1.5 kpc. The density was calculated 
by assuming pressure equilibrium between the two temperature compo¬ 
nents. The radial range for the Bondi radius is shown by the vertical dashed 
lines. 

here (Fig. □•Fig. a also shows the resulting density profile if the 
nuclear PSF flux is increased by 20%. The peak density drops to 
0.49 ± 0.06 cm^^ and the gradient of the profile inside 0.3 kpc ra¬ 
dius is shallower. Fig. [5] shows two possible lines of best-fit to the 
density profile at the centre of M87. The best-fit power-law slope of 
— 1.2 ±0.2 assumed that the PSF subtraction is sufficiently accurate 
and excluded all points beyond a radius of 0.3 kpc. The shallower 
best-fit slope of — 1.0 ± 0.2 was determined using the increased PSF 
subtraction. The density profile is therefore consistent with pocr^* 
within the Bondi radius. 


3.3 Cooling time and dynamical time profiles 

Multi-temperature structure is expected to form in the ICM when 
the radiative cooling time of the gas becomes short with re¬ 
spect to the gas dynamical time. Within 2 kpc, a second temper¬ 
ature component is significantly detected in the hot gas and this 
is also coincident with the extent of the brightest Ha filaments 
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Figure 5. Deprojected density profile. The nuclear flux contribution has 
been subtracted off using a ChaRT simulation (solid circles). The open cir¬ 
cles show the deprojected density profile if the nuclear flux was increased by 
20%. The radial range for the Bondi radius is shown by the vertical dashed 
lines. 



Radius (kpc) 


jYoung et alJl2002l : ISoarks et alJl2004h . We have therefore calcu¬ 
lated the radiative cooling time and dynamical time profiles to 
determine if this observed structure is consistent with the devel¬ 
opment of local thermal instabilities. The radiative cooling time, 
fcooi = (3/2) nkT/rdh, was determined from the single component 
temperature and density profiles. For the two component spectral 
model, the cooling time is approximately a factor of two longer 
for the hotter component and a factor of 2 — 3 shorter for the 
colder component. The d ynamical time was calculat e d usin g the 
mas s profile for M87 oflRomanowskv & KochanekI ( 1200 ll) (see 
also lChurazov et al]|2008h . Fig.|6]shows the cooling time, the free- 
fall time and the fcool/fff profiles for M87. The radiative cooling 
time drops to 3 ± 2 x 10^ yr at 0.2 kpc and is even shorter for the 
cool X-ray temperature component at (1.1 ±0.8) x 10^ yr. The ra¬ 
tio reaches a minimum of fcool/fff = 16 ± 2 at a radius of 0.9 kpc 
and is comparably low to a radius of ~ 4 kpc. 


3.4 Bondi radius and accretion rate 

The SMB FI should accrete from the surrounding hot atmosphere at 
a rate determined by its mass and the gas density and temperature 
at the radius where the SMBH’s gravitational influence dominates 
llBondilll952h . Under the assumptions of spherical symmetry and 
negligible angular momentum, the Bondi accretion rate is given by 


Mb 

Mgyr-i 


0.012 


kBT \ 
keV ) 



Mbh 

lO^Mp 


2 

, ( 1 ) 


where T the gas temperature, iig is the gas electron density and 
Mbh is the SMB FI mass. An adiabatic index y = 5/3 is assumed. 
The Bondi radius can also be expressed as 


Figure 6. Upper panel: the radiative cooling time profile for a single tem¬ 
perature component model and the dynamical free-fall time profile. Lower 
panel: the ratio of cooling time to free-fall time. 


= 0.031 

kpc \keVy 


( 2 ) 


Recent gas-dynamical and stellar-dynamical black hole mass es¬ 
timates for M 87 differ by a factor of two from Mbh = 6.6 ± 
0.4 X 10^Mw dGebhardt et alj|201ll) to Mbh = 3.5+°;^ x 10® Mq 
dWalsh et al.ll2013l) . Therefore, our estimates of the Bondi radius 
and accretion rate are expressed as a range covering these two val¬ 
ues of the black hole mass. 

For a single temperature spectral model, the temperature 
falls to 0.9l/jyjj keV in the centre of M87 and therefore the 
Bondi radius tb = 0.12 — 0.22 kpc (1.5 — 2.8 arcsec). The elec¬ 
tron density at the Bondi radius ranges between 0.31 ±0.04 cm^^ 
and 0.62 ± 0.05 cm^^, which gives Bondi accretion rates Mb = 
0.1 — O.2M0yr^*. Assuming a standard ~ 10% efficiency and 
accretion onto the SMBH at the Bondi rate, the Bondi accretion 
power is Pb = 0.5 — lx 10“^^ erg s^*. Bondi accretion therefore 
would be sufficient to fuel the mech anical power output from the 
jet of Pjaf = 8®~J X 10 ’ ^^ ergs~^ (eg. [Bicknell & Begdrri^ 19991 : 
IPi Matteo et al.ll2003l : iRaffertv et aDl2006l : iRussell et^ 2013h . If 
we consider only the higher, volume-filling temperature compo¬ 
nent of the two temperature spectral model, the Bondi accretion 
rate would be a factor of ~ 5 lower but still sufficient to power the 
jet. For a nuclear bolometric luminosity LboI ~ 2 x 10“^* erg s^*, 
the Bondi accr etion power is also far greater than the nucl ear ra¬ 
diation losses dRevnolds et al.|[l9^ : IPi Matteo et al.ll2003h . Both 
the accretion flow and the jets are therefore radiatively inefficient 
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dpi Matteo et alj|2003h . although the accretion rate may be much 
lower at the event horizon. 


4 DISCUSSION 

The gas density and temperature profiles within the accretion ra¬ 
dius may in principle reflect the gravitational influence of the cen¬ 
tral SMBH and permit us to distinguish between different accretion 
flow models. M87 is one of only a few systems with a Bondi radius 
that can be resolved by Chandra but the analysis is complicated by 
the bright central point source and jet emission. By stacking short 
frame-time observations of M87 taken after the HST-1 flaring pe¬ 
riod and using a ChaRT/MARX simulation to subtract the PSF, we 
determined the gas density to within the Bondi radius and traced 
the multi-temperature structure on these spatial scales. 


4.1 Hot gas accretion 


Within 2 kpc, two temperature components are significantly de¬ 
tected, which are consistent with constant values of 2 keV and 
0.9 keV down to a radius of 0.15 kpc. No evidence was found 
for an increase in gas temperature within ~ 0.25 kpc radius 
where the SMBH is expected to gravitationally heat the ambi¬ 
ent gas, although this may occur at small radii within 0.1 5 kpc 
(eg. lBrighenti & MathewsriT999l : lOuataert & Naravanllioooll . Ev¬ 
idence for central temperature increases o n this scale has been 
found in several nearby elliptical galaxies (|Pellegrini_eLalJl200^; 


iMachacek et al.l200^ : lHumDhrev et alj2008l : Pellegrini et alj201 

and, in NGC3115, the hotter, volume-filling component of a two 
te mperature model ap pears to increase to ~ 1 keV towards the cen¬ 
tre iwong et al.ll201^ . 


However, the hot atmosphere at the centre of M87 has been 
significantly disturbed by the AGN. Several X-ray cavities are co¬ 
incident with radio jet emission in the core of M87, including one 
at a radiu s of only 0.5 kpc to the SE of the nucleus in th e counter-jet 
direction dYoung et al.l2002l : |Forman et al.l2005Ll2007h . This activ¬ 
ity will generate shocks, stir up the gas and cause departures from 
hydro static equilibrium (eg. iForman et al.l l2005l : IChurazov et al.l 
l2008t) . The counter-jet cavity corresponds to a region of apparent 
low pressure and h igh projected temperature dMillion et alJl201U : 
IWemer et al.ll20l3) . The hotter ~ 2 keV component of the multi¬ 
temperature ICM inside a radius of ~ 2 kpc may correspond to gas 
heated by this recent jet outburst. The multi-temperature structure 
could therefore be due to hotter, shock-heated regions and dense 
clumps of cooling gas but this cannot be distinguished from a more 
homogeneous multi-temperature distribution by the azimuthally- 
averaged temperature profiles. Finer spatial binning of the temper¬ 
ature profile extending in to 1 arcsec radius, excluding the region 
affected by the AGN activity and tracing any azimuthal variation 
may be possible with additional observations if the nuclear flux 
continues to decline. However, with the existing data, the central 
temperature gradient is likely affected by the AGN activity and may 
not reflect the underlying gravitational potential. 

The density profile was extracted using finer spatial bins, 
which extend in to 1 arcsec, and is also more robust to uncertain¬ 
ties in the PSF subtraction. The profile flattens within 3 kpc 
with a ridge of increased density at 1 kpc (~ 12.5 arcsec) cor¬ 
responding to the bright rim of the inner bubble. From the den¬ 
sity and temperature profiles, we calculated the Bondi accretion 
rate Mg =0.1 — O. 2 M 0 yr^* across the Bondi radius located at 
0.12 — 0.22 kpc (1.5 — 2.8 arcsec), depending on the black hole 


mass. For accretion onto the SMBH at the Bondi rate and as¬ 
suming 10% efficiency, the accretion rate is roughly two orders 
of magnitude greater than that required by the mechanical power 
output of the jet. The density profile for a Bondi flow is given 
by pocr^P, where (5 = 3/2 for r<^r^. Close to the Bondi radius, 
on the scales probed by the Ch andra observations, a Bondi flow 
has a shallower gradient, (5 < 1 ( lBondilll9^ : lNaravan & Yil|l994 
lOuataert & Naravanll200Cl) . Excluding the region affected by the 
cavity, the density profile inside the Bondi radius in M87 is con¬ 
sistent with pocr^', which may indicate the onset of a Bondi flow. 
However, the ADIOS and CDAF models are also supported by a 
shallow density profile and cannot be ruled out by the existing ob¬ 
servations. Alternatively, the density profile could reflect the radio 
jet activity rather than the underlying gravitational potential. Chan¬ 
dra observations of NGC4889, one of two BCGs in the Coma clus - 
ter with a black hole mass of 2 x 10*° Mq jMcConnell et al.l201 ih . 
show a central cavity inside 0.6 kpc radius, whic h is comparable to 
the recent outburst in M87 dSanders et al.ll201^ . The density dis¬ 
tribution is consistent with no X-ray emission within this radius 
suggesting that much of the material has been driven out by the jet. 
Although we cannot rule out the different models, it is not clear that 
the hot gas structure across the Bondi radius is dictated by the grav¬ 
itational potential of the SMBH and the classical Bondi accretion 
rate may not reflect the accretion rate onto the SMBH. 


4.2 Hot outflow? 


Faraday rotation measurements can provide constraints on the 
accretion rate much closer to the SMBH on scale s of tens of 
Schwarzschild radii {r^). For M87, 1 Kuo et al] ( |2014|) determined 
an upper limit on the accretion rate at 21rs of 9.2 x 10^'* Mq yr^*. 
This limit is at least two orders of magnitude below the Bondi rate 
and it is consistent with the density profile pocr^*, determined by 
Chandra, extending in to this radius. However, this analysis as¬ 
sumes that the observed Faraday rotation originates in the accretion 
flow, the density profile follows a power-law with j6 < 3/2 and an 
equipartition strength magnetic field. 

Observations of Sgr A* and NGC3115 are consistent with the 
onset of an accretion flow inside the Bondi radius and a reduced 
accretion rate ont o the SMBH. Usin g a 1 Ms Chandra observa¬ 
tion of NGC3115, IWong et ^ (1201 4l) found a shallow density pro¬ 
file, pocr^*, and a hot gas component with increasing temperature 
within the Bondi radi us reflecting the gr avitational potential of the 
SMBH. For Sgr A*, IWang et ^ 1 I 2 OI 3 I) showed that the density 
profile is consistent with pocr^*/^, assuming a temperature pro¬ 
file Tocr^*. Faraday rotation measurements of Sgr A* provide an 
upper limit on the accretion rate at r<100rs 1 Bower et aU l2003l : 


lMarrone~et aDl200fiL I 2 OO 7 I: iMacauart et aPbOO^ When combined 

with the X-ray measurements at the Bondi radius, this suggests that 
the flat density profile covers a broad radial range down to r ~ lOOrj 
indicati ng little material captured at the B ondi radius reaches the 
SMBH dYuan et al.ll2003l:IWang et'Zll2013l) . 

These observations are consistent with numerical simula¬ 
tions and analytical models of hot accretion flows, with accre¬ 
tion rate decreasing with decreasing radius so that the majority of 
the material is lost between the Bondi radius and the event hori- 


zon feg. Ilgumenshchev & Abramowiczl 

I 999 I; IStone et all 

I 999 I; 

IStone & Pringle 

200 ll; IHawlev & Balbus 

2 OO 2 I: lYuan & Bui 

20101; 

[Begelmanll2012 

). The radial density profile becomes correspond- 


( l2012h covering a greatly expanded radial range found density 
profiles described by for r > lOrj, favouring the 
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ADI OS model where the majority of the mass is blown out in a 
wind dYuan et al.ll2012h . 

If the accretion rate onto the SMBH in M87 is similarly re- 
duced from the Bondi rate to Mbh = 9.2 x 10^^ M© yr^* (at 21rs, 
IKuo et aLll2014ll . the accretion power is Pbh = 5 x 10^^ erg s^* , 
where an efficiency of 10% has been assumed. This is two or¬ 
ders of magnitude below the Bondi accretion power and close to 
the mechanical power output generated by the jet /jet — 8(^(3 X 
10^^ erg s^*. The bulk of the mass flow will presumably be lost 


in a wind, which will de posit mass and momentum in the sur¬ 
round ing ambient gas (eg. lProga & Begelmarj|20o3 : l0striker et al.l 
l2010h . This will reduce the inflow, controlling the growth of the 
SMBH and making the accretion flow time variable. Calculation 
of the Bondi rate also assumes that the hot gas at the Bondi radius 
is falling toward the SMBH. Instead, the temperature profile indi¬ 
cates the hot gas dynamics may not be dictated by the SMBH and 
the metallicity structure and ionized gas dynamics suggest outflow. 
The alignment of cool, metal-rich gas with the radio jet and cavity 
axes is consist ent with a jet-driven outflow of gas uplifted from the 
cluster centre llBelsole et alj|2001 ; Matsushita et all 20021 : iMolendil 
I 2 OO 2 I : ISimionescu et al.l 20081 : Kirkpatrick et alj 2011 ). The most 
recent AGN outburst generated the inner radio bubbles 10^ yr 
ago and uplifted rough ly 5 x 10^ Mq of gas to ~ 20 kpc radius 
jSimionescu et alj20()8l) . 

The state of the hot gas in M87 differs from the cold gas 
disk, which is clearly dictated by the SM BH taarms etaPI 19941 : 
iMacchetto et alJ[l9^ : IWalsh et al.ll2013h . The central ionized gas 
disk rotates at velocities from —500 to -f500 k m s~' with a mass 
of 4 ± 1 X 10^ Mp) inside 1 arcsec ra dius jpord et al.l 1 19941 : 
iHarms et'n]|l994 IMacchetto et al.lll997h . However, optical and 
UV observations of extended dusty gas filaments found broad emis¬ 
sion lines with both blue and redshifted non-Keplerian compo- 
nents at several hundred km/s with respect to th e systemic velocity 
jTsvetanov et alJl9^ : lFord & Tsvetanovll 19991) . Several filaments 
appear to connect to the disk and the kinematics of the inner fila- 
ment s are consistent with a bi-directional wind b lowing from the 
disk iSparks et alJll993l : lFord & Tsvetanovlll999h . The preponder¬ 
ance of evidence supports an outflow from the hot gas and the disk 
of cold gas. 


4.3 Cold inflow? 


Most of the circumnuclear gas in the rotating disk is presum¬ 
ably in the form of molecular clouds but currently only tenta¬ 
tive detections an d upper limits typically of several xlO^Mp) 


tive detectrons an a upper Irmits typically ot several xiu Mp) 
have been found teraine & Wikliriq 1 199^ : ICombes et al.l 120071 : 
iTan et alj2008l : [^alome & Combesl2008l) . Loss of angular momen¬ 
tum through collisions between clouds and the wind blowing off 
the disk could allow this gas to flow inwards and fuel the SMBH. 
In the ‘cold feedback’ model, these circumnuclear gas disks can be 
fed by overdense, thermally unstable blobs of gas that cool rapidly 
from the surrounding hot atmosphere and lose angular momentum 
through drag in the ICM to accrete onto the di sk dPizzolato & Sokeil 
l2005l :l ^ulsenll 198^ : IPizzolato & Sokejl2010h . Short central radia¬ 
tive cooling times in the hot ICM are strongly correlated with the 
presence of extended cool gas filaments sugges t ing that these orig¬ 
inate in cooling o f the hot gas l lHeckmarJll98ll : ICowie et alJll983l : 
IPeres et alj 19981) . 

Fig. shows that the hot gas in M87 is multi-temperature 
within 2 kpc, which coincides with the extent of the brightest 
Ha filaments dYoung et ^l2002l : ISparks et al.ll2004l : I Werner et al.l 


120101) . Numerical simulations of thermal instability in the ICM 
suggest that this multi-temperature structure will form when the 
gas cooling time drops to a fac tor of ten greater than the dynami- 
cal free-fall tim e, L,nn| /tfF<10 dMcCourt et alJl2oT2l : ISharma et al.l 
l2012l : IVoit et al]l201^ . The ratio of the radiative cooling time to 
the free-fall time in M87 drops to fcool/hf = 16 ±2 at 0.9 kpc and 
is comparably low to a radius of ~ 4 kpc. This appears inconsistent 
with the c r iterion for multi-temperature gas where fcool/hf^lO but 
IVoit et all d2014l) suggest that AGN feedback can cause the min¬ 
imum ratio va lue to fluctuate and a bro ader range of values may 
be acceptable. IPizzolato & Sokej d2005h show that, in addition to 
short cooling times and a suppressed rate of thermal conduction, the 
ambient density profile must be shallow for cold feedback to occur. 
Overdense, rapidly cooling gas will be prevented from reaching an 
equilibrium position in the ICM and can cool to low temperatures. 
The most significant flattening of the density profile in M87 oc¬ 
curs at ~ 3 kpc radius, which also encompasses the brightest Ha 
emission and multi-temperature ICM structure. 

The radiative cooling time drops to only 10^ yr for the cooler 
X-ray temperature component inside 4 arcsec radius. There are 
clear associati ons between the sof t X-ray emission and the optical 
line filaments dSparks et al.ll2004ll and, whilst the innermost fila¬ 
ments are outfl owing, the outer fila ments appear to be falling in to¬ 
wards the disk dSnarks et al.ll2004l) . These overdense, rapidly cool¬ 
ing gas clouds ca n decouple from the ICM and nearly free-fall to 
the cluster centre dPizzolato & Sokeill2010l) . At 0.2 kpc radius, the 
cool X-ray gas component represents ~ 20% of the total X-ray gas 
mass and could lose angular momentum to the hot gas component if 
this is not rotating. The cooler X-ray temperature component could 
therefore provide a source of cold gas clouds within the Bondi ra¬ 
dius with a mass deposition rate of ~ 0.1 Mq yr^*. This would pro¬ 
duce sufficient fuel to power the most recent radio jet outburst. The 
X-ray mass deposition rate measured by the XMM-Newlon RGS in 
al.l arcmin wide aperture (5.2 kpc) is M = O. 9 M 0 yr^* for gas 
cooling radiativ ely from 1 keV to 0 .5 keV and M < 0.06Mq yr^* 
below 0.5 keV dWemer et al.ll201oh . However, if the coolest X-ray 
gas mixes with the cold Ha-emitting filaments, which may be pro¬ 
moted by the AGN activity, then there could be significant non- 
radiative cooling of the gas dFabian et alj|2002l : ISoker et ^|2004 
ISanders et al.ll201(]|) . Taken together, the data are consistent with 
cold feedback where gas cooling from the hot ICM feeds the cir¬ 
cumnuclear gas disk. 


5 CONCLUSIONS 

Using recent short frame-time Chandra observations of M87 taken 
after the flaring period of the HST-1 jet knot, we have subtracted 
the nuclear PSF to determine the temperature and density structure 
of the hot gas within the Bondi accretion radius. Within 2 kpc, we 
detect two significant temperature components, which are consis¬ 
tent with constant values of 2 keV and 0.9 keV down to a radius of 
0.15 kpc. No evidence was found for the expected increase in the 
gas temperature within a radius of ~ 0.25 kpc, where the SMBH 
is expected to gravitationally heat the surrounding gas. This may 
occur at smaller radii within 0.15 kpc. 

From the temperature and density profiles, we calculated the 
Bondi accretion rate of 0.1 — O. 2 M 0 yr^* across the Bondi radius 
at 0.12 — 0.22 kpc (1.5 — 2.8 arcsec), depending on the black hole 
mass. The density profile flattens inside 3 kpc and then steepens 
again within 0.3 kpc (4 arcsec) to peak at 0.5 — 0.6 cm^^ . Inside the 
Bondi radius, the density profile is consistent with pocr^*. Whilst a 
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temperature rise within 0.15 kpc cannot be ruled out, the absence of 
the expected temperature increase inside the Bondi radius indicates 
that the hot gas dynamics are not dictated by the SMBH and, to¬ 
gether with the shallow density gradient, suggests that the classical 
Bondi rate may not reflect the SMBH’s accretion rate. If this shal¬ 
low density gradient extends inwards to the SMBH, consistent with 
the Faraday rotation measurements for M87, the accretion rate onto 
the SMBH could be two orders of magnitude below the Bondi rate. 
Whilst this is still sufficient to power the recent radio jet outburst, 
assuming an efficiency of 10%, the accretion flow will be signif¬ 
icantly altered if a wind blows out the majority of the mass. The 
uplift of metal-rich gas along the radio jet axis and broad optical 
and UV emission lines in the ionized gas disk support the existence 
of an outflow in the hot gas and the cold gas disk. 

Clear associations exist between the coolest X-ray gas and 
the optical line-emitting filaments at the centre of M87 and sev¬ 
eral outer filaments appear to be infalling towards the disk. Obser¬ 
vations show that the ICM has multiple temperature components 
inside a radius of 2 kpc, the radiative cooling time of the coolest 
X-ray gas drops to 10^ yr here and this rapidly cooling gas is co¬ 
incident with the extent of the bright Ha emission. The density 
profile also flattens significantly at this radius, which is consistent 
with cold feedback models that suggest a shallow density profile 
will prevent overdense, cooling blobs from reaching an equilibrium 
position. The cooler X-ray temperature component could provide a 
source of cold gas clouds within the Bondi radius and sufficient fuel 
to power the most recent radio jet outburst. 
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APPENDIX: SIMULATING THE NUCLEAR PSE 
AI Pileup 

Two or more photons landing on the same ACIS detector region 
within the same frame integration ti me will be d etected as a sin¬ 
gle event. Known as pileup (see eg. iDavi sl l200lll . the photon en¬ 
ergies sum to make a detected event of higher energy. The source 
spectrum hardens and the altered shape of the charge cloud distri¬ 
bution can cause grade migration, where the event is misidentified 
as a cosmic ray and excluded from analysis. Most of the intensity 
of piled up sources should still be recovered. However, when the 
jet knot HST-1 was near its peak, strong pileup produced events 
that exceeded the energy filter cut off and were not telemetered 
to the ground resulting in lost flux dHarris et alj|2009ll . The prox¬ 
imity of HST-1 to the bright nucleus also produced the ‘eat-thy- 
neighbour’ effect where photons arriving within the same frame 


and 3x3 pixel grid are considered part of the same event. The pho¬ 
ton w ith the most energy then determines the recorded location for 
both dHarris et al.ll2009l) . It was therefore important to determine 
the extent of any pileup in the observations used and whether it 
was possible to accurately subtract the nuclear PSF from the sur¬ 
rounding extended emission. 

The pileup fraction was estimated using the ratio of ‘good’ 
event grades to ‘bad’ grades in a region of 1 arcsec radius centred 
on the emission peak. Most of the nuclear emission comes from the 
0.5 — 2 keV energy band which, if piled up, increases the number 
of higher energy events above 2 keV with ‘bad’ grades. We there¬ 
fore used the ratio of ‘bad’ to ‘good’ grades in the 2 — 5 keV energy 
band as an estimate of the pileup fraction but note that this measure 
is dependent on the nuclear spectrum and affected by onboard re¬ 
jection. This analysis was confirmed in section A3 with a MARX 
simulatio n of the nucl ear PSF that included a statistical treatment 
of pileup ( lDavisl200lh . 

Apart from obs. ID 352, all of the observations used have a 
short frame integration time of 0.4 s to reduce pileup. We estimate 
that obs. ID 1808 has the lowest pileup fraction at 6% and the 
more recent observations taken since 2010 have a slightly higher 
pileup fraction of 10 — 15% because the nucleus was a factor of 
1.5 — 2 brighter. Obs. ID 352 has a standard frame integration 
time of 3.2 s and is therefore strongly piled up by ~ 80%. Following 
[Harris et ^ ( l2006h . the brightness of HST-1 was determined using 
a detector-based measure of intensity (keV s^*), which sums the 
photon energies from the event 1 file with no grade filtering. The 
HST-1 brightness drops from 0.3 keV s^* in 2010 to 0.2 keV s^’ 
at the end of 2012. These measurements are all more t han an order 
of mag nitude below the 4keV s^* limit above which [Harris et al.] 
l l2009h suggest that more complex pileup problems occur. The mea¬ 
surements also confirm that HST-1 has returned to the 2000 bright¬ 
ness level before the major flare and pileup is now only mild for 
observations with a frame integration time of 0.4 s. 


A2 Nuclear spectrum 

For each observation, we extracted a nuclear spectrum using a re¬ 
gion of 1 arcsec radius centred on the nuclear emission peak. We 
assume that all emission within 1 arcsec is from the AGN. How¬ 
ever the gas emission is likely to be ^ 5% of the total, depending 
on the slope of the gas density profile. A corresponding background 
spectrum, including the cluster emission, was extracted using a sur¬ 
rounding annulus from 2 — 4 arcsec. The spectra were restricted to 
the 0.5 — 7 keV energy band and appropriate responses and ancil¬ 
lary responses were gener ated. The spect ra were fitted individually 
in XSPEC version 12.8.2 (lArnaud l l996l) with an absorbed power 
law model. All parameters were left free and the modified version 
of the C-st atistic available in XSPEC was us ed to determine the best- 
fit values (ICashlll9'^ : IWachter et al.lll979h . Comparing the obser¬ 
vation in 2000 with those in 2010 to 2013, there has been modest 
variation in the best-fit photon index and the absorption parameter, 
uh (Table [TJ, which was consistently abo ve the Galactic value of 
0.019 X 10^^ cm^^ tealberla et alj|200^ . The nuclear flux varies 
by up t o a factor of a few, whic h is consistent with previous stud¬ 
ies (eg. [Harris et alJl997[ . l2009h . The spectra for obs. IDs 11513 to 
14974 were also fitted simultaneously in XSPEC to determine the 
average best-fit model, weighted by the exposure, for the ChaRT 
simulation. The best-fit parameters for this ‘stacked’ spectrum are 
given in Table (Tj 
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Figure 7. Left: Background-subtracted surface brightness profiles in the energy band 0.5 — 7.0 keV for obs. ID 352, obs. ID 1808 and summed obs. IDs 11513 
to 14974. The open symbols show the predicted contribution to the total emission by the nuclear PSF from the ChaRT simulations for obs. IDs 11513 to 14974 
(red) and obs. ID 1808 (blue). The radial range for the Bondi radius is shown by the vertical dashed lines. Right: Zoom in showing the region of overlap 
between the subarray and full frame datasets. 


A3 ChaRT simulations 

From the nuclear spectral model and the s ource position in m irror 
spherical coordinates, the ChaRT program dCarter et al.ll2003l) and 
the MARX softwar^H were used to simulate the nuclear PSF. ChaRT 
provides an interface to the S AOTrace raytrace code, which is based 
on the most accurate mirror mode l and incorporates the details of 
the support structures and baffles ll.Terius et alJll99.'5h . ChaRT sim¬ 
ulations of the PSF core and i nner wings out to 10 arcsec match 
well with on-axis observations (|jeriusl2002h . The M87 nucleus was 
on-axis in each observation analysed and the nuclear emission is 
only significant compared to the cluster emission to 3 arcsec ra¬ 
dius. Therefore the ChaRT simulations should provide a good re¬ 
construction of the PSF. Multiple ray-traces were used and analysed 
simultaneously to ensure sufficient photons for a detailed compar¬ 
ison with the observations. MARX version 4.5 was used to project 
the ChaRT raytraces onto the ACIS-S detector and account for the 
effects of the detector response. Although some differences exist 
between the MARX ACIS effective area and that observed, the ex¬ 
cellent agreement between the subtracted profiles in Fig. [T] shows 
that this is not a significant effect. MARX also includes the standard 
effects of the dither motion and residual blur from aspect recon¬ 
struction errors using the DitherBlur parameter. 

From the simulated events file, an exposure-corrected image in 
the 0.5 — 7 keV energy band was generated with spectral weighting 
determined by the nuclear emission model in section A2. A surface 
brightness profile for the simulated PSF was extracted from this im¬ 
age using a series of concentric annuli of 1 arcsec width. A surface 
brightness profile was also generated for each of the M87 obser¬ 
vations from an exposure-corrected image weighted by a spectral 
model for the cluster emission in the core (section |2^ . The sur- 


' See http://space.mit.edu/CXC/MARX/ 


face brightness profiles for the obs. IDs 11513 to 14974 were then 
stacked to produce an average profile. The short frame integration 
times of these observations, and obs. ID 1808, required the use of 
a subarray to reduce the readout time and avoid losing data. The 
resulting surface brightness profiles only extended to a radius of 
15 arcsec and therefore, despite strong pileup in the centre, obs. ID 
352 was used to trace the cluster gas to larger radii. 

The surface brightness profiles for the M87 observations and 
the PSF simulations are shown in Fig. |3 The nuclear emission 
clearly dominates over the extended cluster emission inside a ra¬ 
dius of 2 arcsec. The profiles appear consistent beyond a radius 
of ~ 5 arcsec as expected if all variable point sources have been 
excluded from the cluster emission. Although there has been a sig¬ 
nificant reduction in the effective area at low energies over time 
due to the build up of contaminant ( [Marshall et alj|2004l) . and this 
has accelerated since 2009, there is still reasonable agreement in 
the cluster surface brightness for the region of overlap between the 
datasets (Fig.|7]right). Flowever, within ~ 5 arcsec, variation in the 
nuclear flux between the observations produces significant differ¬ 
ences. The nucleus was brighter for the stacked observations taken 
from 2010 to 2013 compared to obs. ID 1808, which was taken in 
2000. Obs. ID 352 has a 3.2 s frame integration time and is there¬ 
fore strongly piled up to a radius of 3 arcsec, which reduces the flux 
detected in the 0.5 — 7 keV energy band. 

The stacked 11513 — 14974 observations were mildly piled 
up, which reduced the measured nuclear flux from the central 
0—1 arcsec region over that in the PSF wings. The PSF simula¬ 
tion was repeated usin g MARX to include a statistical realisation 
of pileup ( lDavisll200lh . This piled up simulation was used to esti¬ 
mate the effect on the nucleus model spectrum and showed that the 
0.5 — 7 keV flux in the central region drops by a factor of ~ 1.15 
with little effect in the surrounding annuli. The flux of the PSF sim¬ 
ulation was therefore increased to compensate but, while this was 
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Figure 8. Background-subUacted surface brightness profiles in the energy 
band 0.5 — 7.0 keV for the summed and individual obs. IDs used in this 
analysis. The solid line shows the ChaRT simulation of the nuclear flux. 
The radial range for the Bondi radius is shown by the vertical dashed lines. 


a sufficient correction for the normalization, there was also a low 
level hardening of the spectrum. The soft emission from the hot 
gas that is also present at the ~ 5% level within 1 arcsec radius 
could reduce this effect. However, this uncertainty has limited the 
temperature measurements at small radii in this analysis. 

Fig.[^shows the surface brightness profiles for each of the ob¬ 
servations taken from 2010 to 2013. The nuclear flux varies mod¬ 
estly over this period but beyond 2 arcsec radius the profiles are 
consistent. We note that the increased uncertainty in the 5 — 6 arcsec 
radial bin is due to the presence of a cavity in the X-ray emis¬ 
sion and there are fewer counts in this region. Apart from the cen¬ 
tral source, there do not appear to be any variable point sources 
that have not been removed from the analysis. The ChaRT/MARX 
simulations were tested by comparing the PSF-subtracted surface 
brightness profiles for obs. ID 1808 and the stacked 11513 to 14974 
observations as described in section[2 
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